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I.  INTRODUCTION 


A Joint  Atmospheric  Modeling  and  Chemical  Dynamics  Workshop 
sponsored  by  the  Department  of  Transportation's  (DOT)  Climatic  Impact 
Assessment  Program  (CIAP)  was  held  at  the  National  Bureau  of  Standards  in 
Gaithersburg  on  September  12-13  , 1972.  The  workshop  was  organized  by 
Mr.  Reynold  Greenstone  and  Dr.  Robert  L.  Underwood  of  the  DOT/CIAP  office. 
Host  to  the  workshop  was  Dr.  Milton  Scheer,  Chief  of  the  Physical  Chemistry 
Division  of  the  National  Bureau  of  Standards  . 

The  workshop  was  divided  into  two  separate  panels,  a Chemical 
Dynamics  Panel  chaired  by  Dr.  Ralph  Kummler  of  Wayne  State  University  and 
an  Atmospheric  Modeling  Panel  chaired  by  Dr.  Thomas  D.  Taylor  of  the 
Aerospace  Corporation.  Appendix  A is  a complete  list  of  participants  and  then- 
affiliations  . 

The  prime  objectives  of  the  panels  were  to  review  results  obtained  by 
the  panel  members  that  are  pertinent  to  CIAP  objectives;  to  recommend  models 
that  can  and  should  be  relied  upon  in  addressing  the  CIAP  problem;  to  recom- 
mend best  chemistry  data  for  use  by  the  modelers;  and  to  consider  ways  to 
improve  the  interactions  between  CIAP  modelers  and  experimenters. 

The  following  sections  of  this  report  present  (a)  highlights  of  the  work- 
shop overall;  (b)  the  chemical  dynamics  panel  report;  and  (c)  the  atmospheric 
modeling  panel  report.  Notes  for  the  chemical  dynamics  panel  report  were  pro- 
vided by  Dr.  Thomas  Hard  of  the  Transportation  Systems  Center  who  served  as 
rapporteur  assisted  by  Mr.  Jerome  Pressman  of  Pressman  Enterprises.  The 
report  of  the  panel  begins  with  an  overview  by  the  panel  chairman.  Dr.  Kummler. 
Both  notes  and  an  overview  of  the  atmospheric  modeling  panel  were  provided  by 
the  panel  chairman.  Dr.  Taylor. 
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II.  HIGHLIGHTS  OF  THE  JOINT  ATMOSPHERIC  MODELING  AND 
CHEMICAL  DYNAMICS  WORKSHOP 


ATMOSPHERIC  MODELING 

The  Atmospheric  Modeling  Panel  reached  agreement  upon  requirements 
for  a model  of  the  wake  behind  high-altitude  aircraft,  examined  initial  studies 
of  important  aspects  of  the  dispersion  and  transport  region,  and  learned  de- 
tails concerning  two  separate  approaches  to  modeling  global-scale  circulation. 
Also,  newly  initiated  studies  to  investigate  the  sensitivity  of  global-circula- 
tion models  to  changes  in  input  data  were  outlined.  Many  specific  problem 
areas  requiring  further  attention  were  identified;  these  are  discussed  in  the 
panel  report  (Chapter  IV)  . 

CHEMICAL  DYNAMICS 

Among  the  highlights  of  the  chemical  dynamics  panel  were  the  presenta- 
tion of  new  values  for  certain  reaction  rate  constants  , the  definition  of  certain 
problems  that  need  further  attention,  the  announcement  of  the  availability  of 
O2  photodissociation  computation  programs  , and  an  action  to  coordinate  and 
update  results  of  chemical  rate  determinations  among  members  of  the  CIAP 
chemical  dynamics  community. 

New  Values  for  Reaction  Rate  Constants 

Among  those  reporting  on  new  rate  constant  determinations  were  Scheer 
of  N .B.S  . , Davis  of  Maryland,  Stuhl  of  Ford,  and  Simonaitis  of  Pennsylvania 
State  University.  The  reactions  involved  NO2,  NO,  O,  02,  H,  OH,  HO2  , and 
CO. 
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Definition  of  Problems  Needing  Further  Attention 

Among  areas  needing  further  work  are:  further  verification  of  calcula- 
tions of  NO  photolysis  to  confirm  that  NO  + h^  -►  N + O is  an  important  loss 
mechanism  for  NO  in  the  lower  mesosphere;  resolution  of  the  disagreement  by 
a factor  of  two  for  the  rate  constant  for  NO2  + O NO  + O2;  improvement  of 
O (1D)  rate  constants  to  better  than  the  current  factor  of  two;  determination  of 
the  quantum  yield  of  O (1D)  from  O3  + hv  near  3100  A;  determination  of  the  rate 
constant  for  NO  + HO2  NO2  + OH;  determination  of  the  spectral  quantum 
yield  of  HNO3  photolysis;  and  achievement  of  high-resolution  irradiance  spec- 
tra as  a function  of  altitude,  particularly  in  the  210  nm  window 

The  Availability  of  O2  Photodissociation  Computation  Programs 

Both  J.  Park  (Colorado)  and  R.  D.  Hudson  (NASA)  offered  to  make 
available  their  computational  programs  for  O2  photodissociation. 

Coordination  of  Chemical  and  Photochemical  Rate  Data 

An  agreement  was  reached  between  Johnston  (U.C.  Berkeley)  apd  Garvin 
(N  .B  .S  .)  to  institute  a program  to  collect  rate  data  for  all  chemical  and  photo- 
chemical processes  relevant  to  the  CIAP  problem.  Both  Johnston  and  Garvin 
will  guide  this  effort  with  N .B  .S  . assuming  the  secretarial  role.  Individuals 
with  a special  concern  for  reactions  involving  a particular  species  were  asked 
to  serve  as  an  initial  clearinghouse  for  data  on  those  reactions  and  then  to 
forward  their  results  to  N .B  .S  . 
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III.  SUMMARY  PROCEEDINGS  OF  THE  CHEMICAL  DYNAMICS  PANEL* 

(Ralph  Kummler,  Chairman) 


OVERVIEW  AND  SUMMARY 

The  objective  of  the  chemical  dynamics  workshop  was  to  review  the 
chemistry  and  photochemistry  of  the  natural  and  aircraft -perturbed  stratosphere. 
The  discussions  concerned  chemical  models;  the  evaluation  of  chemical  kinetic 
rate  data;  laboratory  measurements  of  reaction  rates  , absorption  coefficients  , 
and  quantum  yields;  and  solar  flux  measurements.  The  program  began  with  re- 
ports of  work  sponsored  by  the  Climatic  Impact  Assessment  Program  (CLAP)  , 
proceeded  to  other  relevant  efforts  , and  concluded  with  a discussion  of  needs  . 

The  following  summary  of  the  recommendations  of  the  chemical  dynamics 
panel  was  presented  by  panel  chairman  R.  Kummler  to  a joint  meeting  with  the 
atmospheric  modeling  workshop: 

1 . The  set  of  chemical  reactions  being  reviewed  by  the 
CIAP  program  for  determination  of  their  stratospheric 
relevance  is  expanding,  in  order  to  avoid  the  omission 
of  any  potentially  critical  processes  . However,  the 
objective  of  the  chemical  modeling  group  remains  the 
reduction  of  the  chemical  reaction  mechanism  to  the 
smallest  set  needed  to  accurately  describe  stratospheric 
chemistry . 


* Dr.  Thomas  Hard  of  the  DOT  Transportation  Systems  Center  prepared  this 
summary  proceedings  of  the  chemical  dynamics  panel  with  the  assistance  of 
Jerome  Pressman,  Ralph  Kummler,  and  Anthony  Broderick.  Illustrations  and 
some  text  were  kindly  supplied  by  the  speakers  . 
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2.  The  disagreement  by  a factor  of  2 in  measurements  of 
the  rate  of  NO2  + O -*  NO  + O2  at  stratospheric  tem- 
peratures needs  to  be  resolved. 

3.  According  to  Brasseur,  the  predissociation  of  NO  may 
obviate  uncertainties  in  the  downward  flux  of  NO  into 
the  stratosphere;  because  of  its  importance  to  the  ozone 
problem,  this  result  should  be  verified. 

4.  The  quantum  yield  of  O (1D)  in  ozone  photolysis  as  a 
function  of  wavelength  near  310  nm , and  the  reaction 
kinetics  of  0(1D)  , are  critical  problems  . 

5.  The  spectral  quantum  yield  of  HNO3  photolysis  should 
be  measured . 

6.  Modelers  need  absolute  irradiance  spectra  as  a function 
of  altitude,  particularly  in  the  210  nm  window.  High- 
resolution  irradiance  spectra  also  are  desirable.  At  the  • 
steep  310  nm  absorption  edge,  flux  data  obtained  with 
wide -pas sband  filters  are  probably  not  helpful. 

7.  Space  and  time  variations  of  ozone  column  and  solar 
irradiance  are  important,  and  should  be  taken  into  ac- 
count during  the  measurement  of  other  photochemical 
variables . 

REVIEW  OF  JET  AND  WAKE  CHEMISTRY— ^ 

The  present  investigation  is  directed  at  the  nature  and  extent  of  chemical  reactions  in 
the  exhaust  jets  and  wakes  of  high-altitude  aircraft.  The  object  is  to  determine  what 
reactions  might  have  an  important  effect  on  the  domposition  of  hot  gases  as  they  issue 
from  the  engine  into  the  surrounding  stratosphere.  Experimentally  measured  exhaust 
compositions  at  the  exit  plane  of  an  engine  will  not  reflect  accurately  the  final  chemical 
composition  of  the  exhaust  in  the  stratosphere  if  reactions  in  the  wake  arc  important. 


The  objective  of  Task  I,  Problem  Definition,  is  to  identify  essential  features  or  con- 
siderations needed  for  subsequent  development  of  a wake  flow  and  chemistry  model  in 
order  to  define  species  concentrations  and  wake  configurations  at  a time  in  the  wake 
history  when  aircraft-induced  motions  have  diminished  and  are  dominated  by  the  natural 
stratospheric  motions.  The  technical  approach  for  the  wake-chemistry  portion  of 


— This  material  is  taken  from  the  presentation  by  L.  B.  Anderson  of  Lockheed, 
Palo  Alto  Research  Laboratory.  It  is  part  of  a preliminary  draft  report  to 
CIAP  . 
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Task  I has  been  to  screen  the  many  possible  chemical  reactions  to  determine  their 
probable  importance  in  the  engine  exhaust  jet,  the  near  wake,  and  the  far  wake  — i.e. , 
to  as  late  as  20  min  after  passage  of  the  aircraft.  To  examine  the  reactions  in  the 
exhaust  jet,  however,  it  has  been  necessary  for  the  calculations  to  cover  areas  even 
farther  upstream  — i.e. , within  the  engine  nozzle.  This  has  been  necessary  to  examine 
the  rate  histories  for  many  of  the  species  before  they  reach  the  nozzle  exit  so  that  their 
concentrations  at  the  exit  plane  can  be  established.  Many  atoms  and  radicals  are  found 
to  be  present,  and  it  cannot  be  reasonably  anticipated  that  experimental  test  data  on 
representative  turbojet  engines  will  be  capable  of  quantitatively  identifying  all  of  the 
important  species.  Thus,  the  engine  nozzle  kinetics  work  is  a necessary  part  of  the 
chemical  reaction  screening  process.  In  this  approach,  simplifying  realtionships 
have  been  sought,  ouch  as  steady  Btatc  and  partial  equilibrium  relationships  among 
minor  species,  to  discover  when  such  relationships  can  be  used  for  malting  reliable 
predictions.  The  reliability  of  the  simplifying  methods  is  ascertained  by  comparison 
of  the  results  with  fully  coupled  computer  solutions  of  the  reacting  flow  problem. 


FIGURE  1 
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The  emphasis  in  final  development  of  the  chemistry  model  will  be  influenced  by  the 
interface  with  and  requirements  of  the  chemical  modeling  program  for  the  post-wake 
dispersion  regime.  It  is  anticipated  that  the  NO^  concentrations  and  total  NO^  will 
be  first  on  the  list  of  priorities.  It  is  anticipated  that  any  significant  conversion  of 
NO  in  the  wake  to  UNO  will  be  of  major  interest.  Formation  of  sulfate  aerosol  in 

X ^ 

the  wake  and  diffusion  and  transport  regime  will  also  be  of  interest.  It  seems  likely 

that  the  diffusion  and  transport  regime  and  subsequent  climatic  analyses  will  be  less 

sensitive  to  tmeertainlies  in  wake  CO  and  hvdrocarbon  levels  than  to  NO  concentration. 

x 

In  light  of  this,  a simplified  wake-chemistry  model  may  serve  nicely.  If,  on  the  other 
hand,  detailed  information  on  minor  wake  species,  exhaust  generated  radicals,  etc.  , 
is  required,  a much  more  general  and  comprehensive  wake  kinetics  model  will  be 
needed. 

The  scope  of  the  chemical  screening  process  is  broad;  i.e.  , thcrmochemical  reactions, 
photochemistry , and  particulate  generating  reactions  are  all  considered.  Because 
oxides  of  nitrogen  are  thought  to  play  a particularly  important  role  in  the  stratosphere, 
their  reactions  are  given  special  attention,  and  analysis  of  the  photochemical  catalytic 
NO^  and  ozone  histroy  following  passage  of  an  aircraft  through  the  stratosphere  is 
carried  out  to  times  well  beyond  the  20-min  "end  to  wake"  time. 

For  studies  of  thermal  reactions  in  the  engine  exhaust  regime,  the  analytical  approach 
involves  setting  up  a chemical  kinetics  model  using  those  elementary  reactions  that 
play  a role  in  changing  the  concentrations  of  engine  emissions  in  the  exhaust  nozzle 
and  the  near  jet  downstream  from  the  engine  exit  plane.  Numerical  methods  are  used 
to  integrate  the  resulting  governing  chemical  rate  equations  up  to  the  point  where  the 
exhaust  flow  mixes  with  the  ambient  atmosphere.  Initial  conditions  for  the  computations 
arc  compatible  with  available  data  on  engine  exhaust  composition. 

In  particular,  reactions  affecting  the  oxidation  of  CO  to  CO^  and  NO  to  NO^.  as  well 
as  conversion  of  NO,,  to  HNO^  , arc  of  interest.  A model  consisting  of  elementary 
forward  and  backward  reactions  pertinent  to  these, exhaust  chemical  systems  has 
been  used  in  this  study.  Reaction  rate  constants  tire  taken  from  recent  evaluations 
available  in  the  literature.  The  fluid-dynamic  structure  of  the  exhaust  jet  region  was 
obtained  from  analytical  predictions  and  experimental  data  available  from  the 
literature. 
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TABLE  1 


CHEMICAL  REACTIONS  IN  THE  COj/NO^/AIR  SYSTEM 


Reaction 

Number 


Reaction 


Reaction  Rato  Variable** 


n 

E 

Ref. 

0 

10,800 

3-15 

0 

0 

3-15 

0 

24,000 

3-16 

0 

83 , 000 

3-16 

0 

1,930 

3-17 

0 

0 

3-63 

0 

75,580 

3-16 

0 

38,340 

3-16 

0 

65,000 

3-16 

0 

47,400 

3-16 

0 

65,000 

3-16 

0 

56,800 

3-16 

0 

1,080 

3-4 

0 

16,800 

3-18 

0 

1,  100 

3-18 

1 

8,900 

3-18 

0 

5,150 

3-18 

0 

96,000 

3-18 

-2 

0 

3-18 

-1 

118,500 

3-17 

0 

-1,000 

3-18 

0 

1.900 

3-18 

0 

0 

3-64 

0 

1,000 

3-18 

0 

700 

3-18 

0 

1,000 

3-18 

0 

0 

3-19 

1 

2 

3 

4 

5 

e 

7 

6 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 
27 


H ♦ N20  » N2  + OH 
N + OH  - NO  + H 

o ♦ n2o  - n2  + o2 


H + NO, 


2 - NO  ♦ OH 
H02  + NO  = OH  + NO, 

O 

NO  + O 

f'+ 

NO  +'p2  - O + NOg 
NO  + Ntp  = N20  + O 
M + N20\-  N2  + O + M 
CO  + OH  1 


N * 02 
■ O + NO 


M 


H + 02  > 
Oil  + OH 


O + H2  - H + OH 


C02  * H 
OH  + O 

- h2o  + O 


H2  + OH 


H20  + H 


m + h2  - H+  H + M 


M + H +,  OH 


HzO  + M 


M + O. 


O + O + M 


M + H + 02  = H02  + M 
H + HOz  = OH  + OH 

- h2o  + o2 

CMi  + 02 
H + H02  - H2  + 02 
H + H02  = H20  + O 
N02  + OH  = HNOg 


OH  + H02 
O ♦ HO„ 


01  x 10 
21  x 10 
.61  x 10 
41  x 10 
25  x 10 
00  x 10 
44  x 10 
10  x 109 
.10  x 10 
05  x 10 
05  x 10 
30  x 10 
20  x 10 
.20  x 10 
30  x 10 
80  x 10 
20  x 10 
20  x 10 
40  x 10 
50  x 10 
SO  x 10 
50  x 10 
20  x 10 
.00  x 10 
50  x 10 
00  x 10 
90  x 10 


14 


•Reaction  rate  constant  k = A T exp  (E/RT). 
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TABLE  2 


CHEMICAL  KINETIC  SITUATION  TWO  METERS  DOWNSTREAM  OF  NOZZLE  EXIT  PLANE 
FORWARD  RATE  FORWARD  REACTION 


REACTION 

CONSTANT 

RATE 

NET  REACTION  RATE 

COMMENTS* 

fVi 

= J VS 

k 

f 

r*  ■ w [Aj)  1 

rf  - kr  * (A.)*} 
) J 

1 

H + NzO 

at  Nj  ♦ OH 

1.01 

E411 

1.57  E-12 

1.57  E-12 

F 

s 

N + OH 

— NO  4 H 

4.21 

E413 

3.21  E-17 

-3  83  E-17 

B 

3 

O + NjO 

r n2  4 °2 

1.14 

E+6 

1.39  E-15 

1.39  E-15 

F- 

4 

N2+N°2 

— NO  4 NjO 

1 3S 

E-5 

1.94  E-28 

-7  48  E-21 

R 

5 

H 4 NOz 

— NO  4 OH 

2.82 

E414 

1.86  E-10 

1.88  E-10 

F 

• 

HOj  + NO 

- OH  4 NOj 

1 25 

E411 

1.51  E-ll 

1.51  E-ll 

F 

T 

o + n2 

— NO  4 N 

6 93 

E-4 

7.69  E-20 

-3  88  E-17 

R 

8 

NO  ♦ O 

- N 4 °2 

6 73 

E+0 

3 30  E-19 

1 78  E-19 

B 

• 

N02  + M 

= NO  4 O 4 M 

1 41 

E41 

2 72  E-20 

-1  75  E-10 

R 

10 

NO  + Oj 

= o 4 no2 

4.20 

E41 

3 93  E-16 

-4.01  E-12 

R 

11 

NO  + NO 

= NjO  4 0 

9 01 

E-4 

5 06  E-23 

-1  94  E-18 

R 

12 

N20  + M 

= N2  ♦ O 4 M 

8 09 

E+l 

2 56  E-18 

-4  42  E-14 

R 

13 

CO  + OH 

= C02  4 H 

3 17 

E411 

2 32  E-8 

1 85  E-8 

B 

14 

H + 02 

— OH  4 O 

3.11 

E410 

3 26  E-7 

-2  30  E-9 

FASTEST  REACTION 
BUT  VERY  NEAR  PE 

IS 

OH  ♦ OH 

= h2o  4 O 

3.53 

E+12 

5.91  E-8 

-8.48  E-12 

PE 

18 

o + h2 

CH40H 

1.57 

E411 

8 06  E-9 

9 74  E-10 

NEAR  PE 

17 

H2  + OH 

:h204H 

1 45 

E+12 

4.65  E-8 

5.63  E-9 

NEAR  PE 

16 

H2  + M 

— H 4 H 4 M 

2.22 

E-8 

3 96  E-24 

-2  26  E-10 

R 

18 

H + OH  4 M 

= H20  4 M 

1 54 

E+l  7 

3 80  E-9 

3 78  E-9 

F 

20 

02  + M 

*=  O 4 o 4 M 

1.85 

E-12 

5 26  E-28 

-3  46  E-12 

R 

21 

M + H + Oj 

= HOj  4 M 

2 54 

E+l  5 

1.91  E-8 

1.91  E-8 

F 

22 

H + H02 

= OH  4 OH  . 

9. 175 

E+13 

1 24  E-8 

1.24  E-8 

F 

23 

OH  + H02 

- H2o  4 o2 

5 90 

E+12 

3 90  E-10 

3.90  E-10 

F 

24 

o + ho2 

= OH  4 Oj 

2.95 

E+13 

3.12  E-9 

3. 12  E-9 

F 

25 

H + H02 

= H„  4 O, 
2 2 

1 73 

E+13 

2 34  E-9 

2.34  E-9 

F 

28 

H + II02 

= h2o  4 O 

5 90 

E+12 

7 98  E-10 

7.98  E-10 

F 

27 

NOz  4 OH 

= hno3 

2.90 

E+12 

1.00  E-12 

-1.82  E-17  \ 

PE 

•F  - FORWARD  REACTION  DOMINATES. 

R - REVERSE  REACTION  DOMINATES. 

PE  - FORWARD  RATE  - REVERSE  RATE,  REACTION  IN  PARTIAL  EQUILIBRIUM 
B - BOTH  FORWARD  AND  BACK  RATES  SIGNIFICANT  AND  NONEQUAL. 
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In  particular,  reactions  affecting  llie  oxidation  of  CO  to  CO^  and  NO  to  NO^.  as  well 
as  conversion  of  NO  to  HNO„  , arc  of  interest.  A model  consisting  of  elementary 

Z o 

forward  and  backward  reactions  pertinent  to  thcset exhaust  chemical  systems  has 
been  used  in  this  study.  Reaction  rale  constants  are  taken  from  recent  evaluations 
available  in  the  literature.  The  fluid-dynamic  structure  of  the  exhaust  jet  region  was 
obtained  from  analytical  predictions  and  experimental  data  available  from  the 
literature. 

The  major  conclusions  reached  in  the  investigation  are  as  follows: 

9 In  an  afterburning  engine,  nonequilibrium  effects  are  an  important  asj>ect  of 
the  chemistry  in  the  engine  exhaust  expansion  nozzle  and  exhaust  jet  regime. 
There  is  substantial  oxidation  of  CO  to  COg  in  the  nozzle  and  jet,  and  the 
concomitant  production  of  H atoms  maintains  the  concentrations  of  H.  O,  and 
> OH  well  above  their  equilibrium  values.  These  nonequilibrium  concentrations 
persist  from  the  exhaust  jet  out  to  the  regions  where  the  jet  mixes  with  the 
ambient  atmosphere,  as  shown  in  Figure  2. 

« Oxidation  of  NO  to  NO  is  not  a significant  factor  in  the  exhaust  regime. 

The  NO^/NO  ratio  is  maintained  well  below  its  equilibrium  value  (~  10  ) 

by  fast  reactions  such  as  the  following: 

N02  + H - NO  + OH 
N02  + M NO  + O + M 
N02  + O — NO  + 02 

• Although  OH  concentrations  are  high  in  the  jet,  the  small  NO,,  concentration 
in  the  exhaust  region  prevents  significant  formation  of  HNO„  from  the 

reaction 

N02  + OH  — HN03 

• In  ground  tests  the  exhaust  gases  are  maintained  at  substantially  higher  tem- 
peratures for  longer  residence  than  at  high  altitude.  As  observed,  the  chemical 
kinetics  model  predicts  the  occurrence  of  significantly  more  CO  oxidation 

in  the  exhaust  jet  for  an  engine  under  static  lest  conditions  at  sea  level, 
than  it  predicts  for  stratospheric  flight. 

® The  subject  of  chemical  kinetics  in  turbulent  flows  is  considered.  To  a large 
extent,  chemistry  is  found  to  be  "diffusion  limited"  in  the  exhaust  jet  core 
and  "recation  rate  limited"  in  the  wake.  The  simple,  "well  mixed"  kinetics 
model,  which  ignores  diffusion  and  turbulence  in  the  exhaust  jet,  appears  to 
agree  reasonably  well  with  ground-test  data  on  CO  profiles.  This  may  mean 
that  the  net  influence  of  turbulence  on  the  chemistry  is  small. 
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Fig.  2 . Centerline  Concentration  Profiles  for  GE-4  at  Maximum  Power 


LOCKHEED  PALO  ALTO  RESEARCH  LABORATORY 
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• In  the  exhaust  jet,  ratios  between  H,  ,0,  and  OH  appear  to  be  nearly  constant, 
independent  of  elapsed  time.  Effects  tending  to  alter  the  ratios,  6uch  as  O- 
atom  removal  by  methane,  appear  to  be  small.  This  allows  6tcady-state 
approximations  to  be  used  for  HOg/OH,  N02/N0,  and  S03/S02  ratios,  which 
are  6 mall . 

• In  the  jet  mixing  region  and  the  near  wake,  the  possibility  of  aerosol  formation 
by  the  following  reaction  pair  cannot  be  ruled  out: 

S02  + O + M — S03  + M 

so3  + h2o-  h2so4 

The  magnitude  of  the  effect  is  expected  to  be  highly  dependent  on  the  specific 
cooling  and  dilution  rate  histories  involved.  All  other  reactions  that  might 
lead  to  particulate  formation  appear  to  be  too  slow  to  contribute  significantly 
in  the  wake  regime. 

• Addition  of  unburned  exhaust  hydrocarbons  to  the  stratosphere  does  not  appear 
to  be  significant.  Catalytic  chain  effects  of  hydrocarbons  on  ozone/NO^  cycles 
degrade  the  hydrocarbons  and  are  unimportant.  Removal  of  NO^  by  hydro- 
carbons in  the  stratosphere  is  an  important  area  for  study,  but  the  contri- 
bution from  engine  exhaust  hydrocarbons  is  dwarfed  by  tho  ambient  methane. 


LOCKHEED  PALO  ALTO  RESEARCH  LABORATORY 
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REVIEW  OF  CHEMICAL  DYNAMICS  FROM  THE  DIFFUSIVE  TRANSPORT 
SCALE  TO  THE  GLOBAL  SCALE  1/ 


I.  ASILOMAR  STATUS 

A.  Data 

. List  of  Chemical  and  Photochemical  Reaction  Processes  - 
Partially  Ordered  for  Importance  and  Directed  Mainly 

at  03. 

. Preliminary  Estimate  of  Uncertainties  in  Rate  Coefficients 
for  Sensitivity  Studies 

, No  Comprehensive,  Up  To  Date,  Compendium  of  Evaluated 
Atmospheric  Reaction  Rate  Constants. 

. Recommendations  for  Additional  Laboratory  l.'ork. 

. Recommendations  for  Updated  Atmospheric  Composition, 

Solar  Flux,  and  Optical  Data. 

B.  Models 

Operate,  to  date,  upon  simplified  forms  of  the  more  general 
set  of  kinetic  equations: 

+ V-V  C.  = F.(C)-Ri(C)Ci  + S.(x,t)  + fLlx.t)  , (1) 

v/here  C..  = C.(x,t)  is  the  mean  concentration  of  snecies  i at  position 
x and  time  t, 

•*  -V 

V = V(x,t)  is  the  velocity 


2/ 

This  material  is  taken  from  the  presentation  by  R.  Gelinas  of  Lawrence 
Livermore  Laboratory  in  which  he  reviewed  the  status  of  CIAP  chemical 
modeling  as  of  the  Asilomar  Conference  held  on  11-14  April  1972  (pro- 
ceedings available  from  the  CIAP  Program  Office) , developments  since 
Asilomar,  and  present  and  future  needs. 
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F..(C)  = local  production  rate  by  chemical  and  photochemical  reactions 

R..(C)C..  h local  destruction  rate  by  chemical  and  photochemical 
reactions 

S^(x,t)  = local  net  source  term 

H.j(x,t)  = terms  representing  departures  from  mean  variables 

Classes  of  Models  Receiving  Majority  of  Attention  at  Asilomar 

3C. 

• Fi(C)-Ri(C)Ci  + S.(t«)  (2) 

(often  set  S^/at  = 0 and/or  S^(x,t)  = 0) 

•+  -V 

. V* V = F- (C)-R^ (C)C^  (one-dimension)  (3) 

SC.  -*  -*■ 

H C.  = S.(x,t)  + Hi(x,t) 

(one,  two,  or  three  dimension)  (4) 

(often  set  S^x.t)  = 0 and/or  l^(x,t)  = 0) 

. Mechanism  Modelling  Onerates  on  Sets  of  Chemical  and  Photochemical 

• — ► -v 

Reactions  to  be  included  in  F^ (C)-Ri (C)C^  Terms. 

, Inhomogeneous  Chemical  Modelling  Places  Emphasis  Upon  H^(x,t)  terms 

C.  Standard  Atmosphere  and  Solar  Fluxes 
. U.S.  Standard  Atmosphere  p,  p,  T vs.  z 

. concentrations  of  minor  species  vs.  z badly  lacking  for  flO,  NO^, 
sulfur  compounds,  and  other  minor  species  which  may  possibly 

o o 

absorb  from  1700  A - 3000  A 

. microscopic  absorption  cross  sections,  o (x,p,  T)  - manv  lacking 

a 

high  resolution,  particularly  in  line  wings 
e.g.  02  - Schumann-Runge  Bands 
NO  - Bands 

o 

. solar  flux  outside  earth's  atmosphere  lacks  resolution  to,  say,  1 A 

. solar  flux  vs.  altitude  - a number  of  routines  internal  to  model 
codes  exist  - apparently  none  incorporate  high  resolution  data  for 
02  Schumann-Runge  bands  of  Ackerman,  et  al.  The  reduction  factors 
Kockarts  are  not  particularly  useful  for  general  modelling  needs. 
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D.  Aerosols 


. Major  concerns  are 

i)  Radiative  properties 

ii)  Surface  chemistry  - effect  on  stratospheric  chemistry 
directly  because  of  capacity  as  sinks 

. Qualitative  considerations  exceed  quantitative  work,  to  date 


II.  DEVELOPMENT  SINCE  ASILOMAR 

A.  Data 

. Chapman  reaction  rates  improved 

. NBS  updating  earlier  compendium  of  evaluated  atmospheric 
rate  constants  (numerous  evaluators  included) 

NBS  considers  only  evaluated  data 

. ILL,  in  cooperation  with  NBS  and  individual  users,  maintainino 
compendium  of  evaluated  data  and  other  data  being  used  bv 
various  modellers  (particularly  in  absence  of  evaluated  rates). 

, LLL  states  what  they  are  using  in  LLL  codes  at  a given  time 

for  radiative  and  chemical  rate  data. 

. Much  of  hierarchy  in  Asilomar  table  3.1  has  been  documented 
with  regard  to  data  sources  and  data  values  with  uncertainties. 

Also  "consistency  flagging"  is  underway  for  sensitivity 
studies. 


B.  Models 


at 


Fi(C)-R.(C)C.  + S(t) 


(5) 


wi th  i ) 
ii) 

iii) 

1v) 


current  evaluated  data 

high  resolution  (in  Schumann-Runge  bands),  time- 
dependent  solar  flux 
time  varying  sources 

completely  general  solution  method  for  stiff  systems 


r3  K-  ^ f 3 

. + v.VC.pJd3x< 


Fi(C)-R.(C)C1 


(6) 


with  i)  local  (x,t)  solar  fluxes 
ii)  current  evaluated  data 
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most  significant  development  is  formulation  of  dynamic,  kinetic 
models  (one  space  dimension) 

f£  + V-V  C.  = F.(C)-R.(C)C.  + Si(x.t)  (7) 

this  development  represents  an  emergence  from  an  era  of  exceedingly 
macroscopic  models  which  reauired  massive  phenomenological  inter- 
pretation for  detailed  (or  local)  information. 

. Likely  directions  for  future  efforts  brought  about  by  dynamic- 
kinetic  modelling  will  be  discussed  in  Section  III. 

C.  Standard  Atmosphere  and  Solar  Flux  / 

. Are  a number  of  routines,  which  can  calculate  solar  flux  for 
arbitrary  time,  altitude,  latitude,  and  longitude;  given  column 
densities  of  minor  constituents  and  standard  atmospheric  profile, 
for  major  consti tuents-particularlv  0^.  At  least  one  routine 
has  high  resolution  of  Schumann-Runge  bands  (from  data  of  Ackerman, 
et  al). 

. Is  ongoing  theoretical  and  experimental  work  to  improve  treatment 
of  1 ine  profiles  in  0^  bands.  This  could  possibly  be  significant 
effect  in  transparency. 

III.  UPDATING  OF  HEEDS 

A.  Data 

. numerous  needs  in  chemical  rate  constants  - specific  reactions 
evident  in  compendia  and  Asilomar  proceedings 

1)  low  temperature  data  needed  for  large  no.  of  reactions 

ii)  more  precise  data  needed  for  many  reactions  - greatest 
priority  attached  to  reactions  most  affecting  odd  oxygen 
and  to  reactions  associated  with  stratospheric  sink 
mechanisms  for  species  affecting  odd  oxygen 

. radiative  data 

i)  need  more  resolution  in  absorption  cross-sections,  a (v) 

cl 

and  especially  in  quantum  yields, 

ii)  need  temperature  and  pressure  dependence  of  oa(v), 
especially  in  bands. 

iii)  need  better  treatment,  theoretically  and  experimentally 
of  line  profiles,  especially  in  wings  of  lines  (already 
mentioned  for  solar  fluxes) 

i v )  only  limited  evaluation,  to  date 
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B.  Models 


. sensitivity  to  space  and  tine  dependent  diffusion  coefficients, 
sources,  and  net  removal  processes  (sinks)  can  be  as  great  as 
sensitivity  to  major  catalytic  cycles  in  dynamic-kinetic  model. 

, the  above  point  can  create  a situation  whore  weak  mechanisms  in 
the  reaction  rate  hierarchy  take  on  significance  of  major  catalytic 
cycles  if  the  weak  mechanism  is  associated  with  a removal  process 
for  something  like  dO  - thus  may  need  to  adjust  our  existing  inter- 
pretation of  a reaction  hierarchy  to  accommodate  coupled  dynamic- 
kinetic  mechanisms. 

. the  above  points  amplify  need  to  understand  chemistry  of  aerosols 
as  sinks  for  important  snecies 

. the  above  points  amplify  the  need  for  more  accurate  determinations 
of  NO,  0^,  and  tiO^  profiles  with  seasonal  variations  as  well  as 

for  better  knowledge  of  specific  source  and  removal  agents. 

• all  of  the  considerations  so  far  advanced  for  ozone  kinetics  and 

dynamics  should  extend  to  cncomnass  the  other  species  and  solar 
spectral  regions  associated  with  earth's  energy  budget. 


CHEMICAL  MODEL  PRESENTATIONS 

The  following  paragraphs  summarize  the  presentations  of  several 
chemical  modelers  at  the  joint  session  of  the  atmospheric  model  and  chemi  • 
dynamics  panels  . 

G.  Brasseur  (Institute  of  Aeronomy,  Brussels) 

Subjects  covered  were  the  aeronomy  of  NO  and  NO2  , including  pi 
dissociation  of  NO.  A line-by-line  calculation  in  the  NO  6-bands,  made 
necessary  by  the  overlying  O2  Schumann -Runge  bands,  seems  to  result  in 
net  source  of  NO  above  the  stratosphere. 

E.  Hesstvedt  (University  of  Oslo,  Norway) 

Hesstvedt described  a model  with  ozone,  water  vapor,  and  nitroge 
oxide  chemistry  coupled  with  vertical  and  meridional  transport.  His  abstrac 
follows: 
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A theoretical  meridional  model  of  the  0^  layer  is  pre- 
sented. Two- dimensional  transport  by  eddies  and  mean  motion  is 
considered  together  with  photochemical  reactions  involving  0,  N, 
and  H.  The  model  is  used  to  evaluate  the  effect  of  increased 
contents  of  H^O  and  of  N0x  in  the  stratosphere.  It  is  found 
that  a doubling  of  stratospheric  humidity  will  reduce  the  total 
amount  of  0,  by  less  than  1% , v/hile  a doubling  of  NO  will 

* J X 

result  in  an  lQ%  reduction  for  middle  latitudes  and  summer.  A 10? 
increase  in  NO  will  reduce  the  total  0,  by  about  2.8?.  The 
relation  between  UV  radiation  and  total  0^  is  described,  in 
particular,  for  wavelengths  of  biological  interest  (290  to  320  nm) . 
A 1?  decrease  in 'total  0^  will  result  in  a 2%  increase  in 
erythc-mogenic  UV  radiation. 


T.  Chang  (Lawrence  Livermore  Laboratory) 

Chang  discussed  the  derivation  of  vertical  transport  from  chemistry, 
and  vice  versa.  In  a one -dimensional  ozone  chemical  model,  he  has  inserted 
the  stratospheric  NOx  injection  recently  estimated  by  Foley  and  Ruderman  for 
the  1961-62  atmospheric  nuclear  weapons  test  series.  He  finds  that  the  effect 
is  spread  over  a period  of  about  three  years  , with  a peak  reduction  of  8%  of  the 
total  ozone  column  occurring  about  a year  after  the  tests.  This  suggests  that 
the  1960-66  ozone  data  require  more  careful  statistical  analysis  to  detect  pos- 
sible effects  on  the  ozone  of  NOx  from  nuclear  weapons  tests  . 

T.  Shimazaki  (NOAA,  Boulder) 

Shimazaki  described  work  on  a two-dimensional  photochemical  model. 
This  model  was  originally  developed  for  the  mesosphere,  and  is  now  being 
modified  for  the  much  longer  characteristic  times  of  stratospheric  chemistry. 

T.  London  (University  of  Colorado) 

London  described  preliminary  work  on  a photochemical  model,  with 
three-dimensional  transport  provided  by  the  NCAR  circulation  model  of  Kasahara, 
Washington,  and  Sasamori. 

H.  Tohnston  (U.C.  Berkeley) 

Johnston  described  his  study  of  the  global  stratospheric  ozone  budget. 
He  showed  that  ozone  is  produced  by  C>2  + hv  -*  O + O , followed  by  O + 02  + M 
-►  O3  + M,  much  faster  than  it  can  be  destroyed  by  O + O3  -+  O2  + O2  or  by 
transport  to  the  troposphere.  Small  amounts  of  NOx  can  account  for  the  deficit, 
and  the  distribution  of  NOx  is  as  important  as  its  amount. 
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EVALUATION  OF  CHEMICAL  KINETICS  RATE  DATA 


D.  Garvin  reviewed  the  achievements  and  plans  of  the  N.B  .S  . evalua- 
tion team,  relative  to  the  list  of  reactions  assembled  at  the  Asilomar  Conference. 
He  also  described  a time-shared  kinetics  data  system  that  could  provide 
modelers  with  the  most  recent  measurements  and  evaluations  of  reaction  rates  . 
The  system  was  demonstrated  during  this  meeting.  It  would  cost  N.B.S.  about 
$1,000  per  year  to  maintain,  and  would  cost  the  user  between  $6  and  $15  per 
hour  of  connect  time.  Garvin's  notes  follow: 


TABLE  3 

CIAP  RATE  DATA  SUMMARY 


Table  3a.  Summary  of  Status  of  Asilomar  List 


Symbol* 

Present  Status 

No.** 

12- 

Some  Recommendation  Available 

31 

-12 

P12 

NBS  Plans  for  Fall  72  and 
73 

15 

? 12 

Other  Reactions  (Estimates  Feasible) 

5 

12 

Major  Problems  Using  Existing  Data 

12 

63 

* Marks  next  to  reaction  number  on  attached  list. 

**  Each  reaction  is  counted  only  once  - in  the  first  category 
in  which  it  appears  . 

Table  3b.  NBS  Data  Sheets 


Total 

Asilomar 

Issued  1972 

28 

(13) 

To  Be  Issued  Fall  72 

14 

(10) 

Plans  for  1973 

13 

( 6) 

55 

(29] 
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TABU.  4.  FRACTIONS  ANH  CONnOPNCK  LIMITS 

tor  TiiLik  kati:  constants 


Reaction 

No. 

Reaction 

Confidence  Statement 

1. 

Bjsic  Clupman  Mechanism 

~*i 

Oj  + he  - O + O 

. o * 20% 

— -i 

O + O,  - O,  + O, 

k well  known 

P 3 

O + 0,  + M - 0,  + M 

k well  known 

-——4  • — 

O,  + hr  - 0('u)  + Oj 

o t 1 0% 

11. 

Destruction  of  Odd-Oxyp.en  by  NO* 

■ — 5 — 

NO  + Oj  - NOj  + Oj 

k t 20% 

6 

NOj  + hf  NO  + O 

o,  J i 20% 

7 

NOj  + 0 — NO  + Oj 

k + 20.  -50% 

HI. 

Removal  of  NO, 

Ps 

OH  + NOj  - HNOj 

k factor  of  2, 

•9 — • 


10- 

• — 12' 
13* 

Pm 


— 23 — 


24 


25 


function  of  pressure 
k factor  of  2 


HNOj  + he  — H + NOj 

o factor  of  100 

~ OH  + NOj 

IV. 

NO  Production 

NjO  + 0('D)  - NO  + NO 

k factor  of  2 

- Nj  + 0, 

k factor  of  2 

°*  * t 

NjO  + h*  - Nj  + Of'D) 

'c,  t 20%  ' 

IV  A. 

Hypothetical  Reaction  Affecting  NO 

°tl 

NjO  + hv  — N + NO 

Important  if  ojj  = .01 
opsJjO  + hv  -*  Nj  + O) 

Nj(A)  + O - NO  + N 
Of'D)  + M — Oi\ Y + M 

0(3p)'  indicates  hot  0(5P) 

M + 0(3FV  + Nj  - NjO*  + M 
NjO*  - NO  + N 

NjO*  indicates  excited  NjO 

- N,  + O 
-*  NjO 

Of’,,)  + Nj  - NjO* 

Slow 

V. 

Higher  Oxides  of  Nitrop.en 

NO,  ♦ Oj  - NOj  + O, 

k.factorof  I04  uncertainty  at 
Stratosphetic  T depending  on 
activation  energy  = 7*2  keal. 

„ _ ^ . v 

NOj  + hv  - NO  + O, 

O,  3 l 20;// at  peak 
$ factor  of  10 

NO,  + hv  - NO,  + O 

o Small  at  peak,  possibly  important 
at  lower  wavelength 

• - .Photochemical  absorption 
fc-  Reaction  rate  constant 

cross  section 
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TMU.i:  4 (Coot) 


1 

V.  Higher  Oxides  of  Nitrogen  (Coni.) 

2ll  — 

27 

28 

• — 2V 

P 20 

NOj  + N03  ••  N2  Oj  * 

Nl°s*  ~ NOJ  N0J 

- N,0, 

NjOj  + hi-  — ' 

NjOs  + O - | 

N2Os  + OH  - | 

HjO,  + HOj  - 

k « • 2 in 

_ 

being  worked  on 

23 — . 

NjOj  + HjO  - 

Slow  in  gas  phase 

VI.  Of 1 1,  > and  03 

3-1 

0('„)  + Oj  - Oj  + Oj 

35 

Ol'D)  + M - 0(3p)  + M 

VII.  Reactions  Affecting  OH,  II02 

?3C- 

OH  + O — 11  + 02 

k factor  of  2 

? 37  — 

H + Oj  + M - H02  + M 

k factor  of  2;  probably  only 

r* 

important  11  reaction 

? 38 — 

llOj  + 0 - OH  + Oj 

k not  known 

• — 39 — 

OH  + 03  - HOj  + 02 

k factor  of  103 

— 40  — 

OH!  + 03  - H + 0,  + 02 

k depends  on  vibrational  T ; 

t indicates  vibrational 

excitation 

41 

OH  + CO  - COj  + H 

k factor  of  2 

— 42 

0('{,)  + ll20  -*  OH  + OH 

k factor  of  2 

43 

0('D)  + Cll4  - Cllj  -1  OH 

k factor  of  2 

— 44 

0('n)  + H2  - OH  + H 

k factor  of  2 

? 45  — 

OH  + HOj  - HjO  + Oj 

k of  2x  I0‘° 

• 

seems  high 

46  — 

OH  + OH  - HjO  + 0 

k factor  of  2 

*—  47  ■ — 

H + Oj  - Oi^  + 02 

k factor  of  2 

—48  — 

HOj  + H02  - ll202  + 02 

k factor  of  2 

49  — 

H2Oj  + In-  — OH  + OH 

0 factor  of  2 

? SO 

Hj02  + O - HOj  + OH 

k factor  of  2 

• — 5! — . 

H2Oj  + Oil  - HjO  + HOj 

k factor  of  2 

52  — 

HOj  + NO  - OH  + NOj 

k probably  high 

53 

HOj  + CO  — OH  + CO, 

k low.  factor  of  10 

54  — 

OH  + CH4  - HjO  + Cll3 

k factor  of  2 

55  — 

OH  + Hj  - HjO  + H 

k factor  of  2 

© 

OH  + h-  - 0 + H 

0 very  uncertain 

• © 

HOj  + he  - H + Oj 

0 very  uncertain 

? 58 

Cllj  t Oj  HI  - CHjOj  + M 

‘<59) 

ClljOj  + NO  - ClljO  + NOj 

(60) 

CHjOj  t ClljOj  - ClIjO  + CHjO  + Oj 

?<m) 

CIIjO  + Oj  - ClIjO  + HOj 

62  — 

CH}0  + Oil  - II CO  + IIjO 

k factor  of  2 

P « 

HCO  + Oj  - CO  + HO, 

k high; uncertain  but 

sufficient 
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TABLE  5 - 

N13S  Chemical  Kinetics  Data  Survey 

I.  Evaluations  of  Rate  Constants  and  Photochemistry 
Published  in  MBS  Reports  10692  and  10828 


A 

0 + 03  - 202 

HM02  + H ->  products 

°2(1Ag)  + M -*■  02  + H 

HN02  + HO  ->  H20  + N02 

02(1Eg+)  + M 02  + M 

HN02  + 0 -»■  HO  + N02 

A 

0^  + h v (u.v. ) 

hno3  + H ->  ho  + hno2 

0^  + hv  (visible) 

HN03  + HO  -*■  H20  + N03  fa 

A 

H + 03  ->  HO  + 02 

HN03  + 0 -*■  HO  + N03 

A 

HO  + 0^  -+  products 

HN03  + hv  A 

A 

ho2  + ho2  + h2o2  + o2 

H + N02  ->  HO  + NO 

A 

H20  + 0(1D)  f 2110 

H20  + NO  + N02  2 HNO 2 

A 

H202  + hv 

Ho0  + No0  - -*■  2 HN0_  A 

2 2 b 3 

A 

NO  + 03  + N02  + 02 

^*2^2  + ^ products 

N02  + 0 + M -*•  NO^  + M 
N02  + O3  N03  + 02 

A 

N20  + 0(1D)  -*■  N + 02 

A 

-»■  2N0 

UNO  + H -►  H2  + NO 
HNO  + HO  H20  + NO 

A = Asilomar  list 
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TABLE  5 (Cont) 


III. 

A 

A 

A 

A 


II.  Evaluations  of  Kinetic:;  and  Photochemical  Data 
to  be  issued,  Fall  1972. 

02  + M ->■  20  + M 

H202  + H + 

A .II202  + HO  -►  HO 2 + U20 
A no2  + 0 NO  + c2 
(,\  02  + h v 

NO  + hv 
A N02  + hv 

A N2°5  + hV 

11N02  + hv 

A OC1!))  + FI, 

m = o2,  k2,  o3>  h2,  ch^,  (n2o,  ji2o) 


Hate  Evaluations  Planned  for  1973 


0 

+ 

02  + M ->  03  + M 

M°2 

+ 

110  3 

NO  + N02  + 0 

0 

+ 

N„Or  

M0„ 

+ 

M0o  + 

M ->-  No0c  + M 

2 9 

2 

3 

0 

+ 

ch2=ch2  -»■ 

■ N03 

+ 

mo3  -> 

M 0 2 + M 0 2 + i 

H?0  + hv  H + HO 


MO  + M03  ->  N02  + N02 


N20  + hv+  II 2 + 0 ( XD ) 


HO  + N02  + M -*■  HM03  + M A 


N + NO 


CHO  + 02  -►  CO  + H02 


A 


CH20  + hv->  CHO  + H 


CO  + HO 2 ->■  C02  + HO 


A = Asilomar  list 


TABLE  5 (Cont) 


IV. 

Other  Possibilities. 
Evaluation 

Evaluate  or  Validate  Existing 

A 

H + 02  + M -»  HO 2 + M 

CH3  + 02+  M.  ->■  ClI30p  + M 

A 

HO  + HO 2 -»•  H20  + 02 

CH2-CH2  + HO  -*■  

A 

HO  + 0 ■>  H + 02 

ch2=ch2  + no2  -v  

>I02  + ‘NO  ->  HO  + N02 

ch2=ch2  + o3  - 

A 

H02  + 0 HO  + 02 

NO  + N02  N203 

N + NO  N2  + 0 

N^Or  + M N0o  + N0_  + M 

d d J 

N + 02  -►  NO  + 0 

N0?  + N02  + M N20i]  + M 

NO  + 0 + M -*•  N02  + M 

N03  + hv 

CHpO  + 0 + CHO  + HO 

V.  Information  on  Rate  Data 

A.  "Selected  Rate  Constants  for  Chemical  Reactions  of 
Interest  in  Atmospheric  Chemistry,"  NBS  Report  10867, 
June  1978 

260  evaluations  or  recommendations  for  rate  constants 
of  180  reactions. 

B.  Possibilities  for  1973 

1.  Semiannual  issuance  of  "Selected  Rate  Constants" 

2.  Time-shared  Rate  Data  System 

- Retrieval  and  rate  constant  calculation 

- "Selected  Rate  Constants" 

- "New  measurements" 
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KINETICS  RATE  MEASUREMENTS 


NBS  Rate  Measurements 

M.  Scheer  reported  the  following  efforts  by  his  group  at  NBS,  who 
also  form  the  nucleus  of  evaluators  in  Garvin's  program.  A.  Bass  and  A. 

Laufer  are  studying  the  ultraviolet  absorption  spectrum  of,  and  equilibrium 
between,  NO2  and  N2O4.  Their  results  at  298  K agree  to  within  10-15%  of 
those  of  Blacet  and  Hall;  they  plan  to  extend  the  temperature  range  down  to 
230  K. 

M . Kurylo  has  studied  H + O2  + M HO2  + M by  resonance  fluorescence 
in  the  range  200-400  K and  10-400  torr.  With  helium  as  third  body,  he  obtains 
k (He)  = 6.7  x 10"33  exp  (238/T)  cm6  molec  ”2  sec-1/  and  the  following  relative 
third -body  efficiencies: 

CH4  N2  He  Ar 

15.7  3.4  1.0  1.0 


His  estimated  error  is  +15%.  He  has  also  measured  the  rate  of  H + N02-"* 

HO  + NO  (5  x 10"11  cm3  molec-1  sec-1  at  298  K)  , and  plans  to  study  the 
reactions  of  ozone  with  H,  O,  and  HO,  all  by  resonance  fluorescence. 

R.  E.  Huie,  working  with  J.  Herron  at  NBS  and  D . Davis  at  the 
University  of  Maryland,  has  measured  the  rate  of  O + O2  + M -»  O3  + M by 
resonance  fluorescence  in  the  range  200-346  K and  50-500  torr.  With  Argon 
as  third  body,  they  obtain  k (Ar)  = 6.6  x 10  "3  3 exp  (511/T)  cm3  molec”2 
sec”1 , estimating  their  error  as  + 15%.  At  218  K , k (N2)  = 1.2  x 10-33,  a 
factor  of  2 lower  than  previous  estimates  . They  have  also  measured  the  rate 
of  O + NO2  "*02+  NO  in  the  range  230-340  K,  obtaining  k = 9 . 1 x 10“12 
cm3  molec"!  sec"1  independent  of  temperature,  a factor  of  2 faster  than  pre- 
vious estimates  . For  the  three-body  reaction  O + NO2  + M + NO3  + M,  they 
obtain  k (Ar,  298  K)  = 2.1  x 10"31  and  k (Ar,  263  K)  = 4 x 10"31,  observing 
a complex  pressure  dependence  below  250  K. 

Herron  and  Huie  are  studying  the  reactions  NO  + O3  + NO2  + O2  and 
NO2  +03"*  NO3  + O2  by  stopped-flow  mass  spectrometry. 

University  of  Maryland  Rate  Measurements 

D.  Davis  (University  of  Maryland)  reported  in  further  detail  on  the 
measurements  of  O + O2  + M , O + N02/  and  O + NO2  + M described  above, 
and  pointed  out  that  these  results  can  change  odd -nitrogen-odd-oxygen  models 
by  a factor  of  4 . For  H + O2  + M + HO2  + M , Wong  and  Davis  have  obtained 
k (Ar)  = 6.75  + .80  x 10”33  exp  [(685  + 128)/RT]  cm3  molec"2  sec“l , with 
relative  third -body  efficiencies: 

Ar  He  H2  N2  CH4 

1.0  0.97  3.0  2.8  21.5. 
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Davis  described  the  controversy  over  the  rate  of  CO  + HO2  CO2  + 
HO,  with  estimates  ranging  from  k=  10“24  (Baldwin  et  al)  to  k = 10~12 
(Westenberg  and  DeHaas)  , and  reported  his  own  measurements.  He  sets  an 
upper  limit  of  10"19  to  this  rate;  obtains  a rate  of  1 0 — 1 5 for  HO2  + SO2 
HO  + SO3;  and  sets  a lower  limit  of  5 x 10 2 to  the  rate  of  HO2  + NO 
HO  + NO2  . 

Pennsylvania  State  University  Rate  Measurements 

R.  Simonaitis  reported  an  upper  limit  of  3 x 10 “19  for  HO2  + CO  and 
a lower  limit  of  5 x 10“13  for  HO2  + NO,  in  agreement  with  Davis  , though 
different  techniques  were  used. 

Ford  Motor  Company  Rate  Measurements 

F.  Stuhl  reported  Ford  Motor  Company's  program  to  measure  the  re- 
actions of  OH  radicals  with  CO  , H2 , D2  , NO , NO2  , NH3  , H2S  , HxCy , and  O3  , 
and  also  the  reactions  forming  OH  by  photolysis  of  H2O,  by  reactions  of  O('D) 
with  hydrogenic  compounds  , by  the  reaction  of  H with  NO2  , and  by  the  reac- 
tions of  HO2  with  CO  and  NO.  In  all  of  these  studies,  OH  formation  and  decay 
are  observed  by  resonance  fluorescence.  The  results  to  date  are  summarized  in 
the  accompanying  abstract. 

ON  THE  KINETICS  OF  Oil  IN  THE  PULSED  VACUUM-UV 
PHOTOLYSIS  OF  MIXTURES  OF  CO,  02 , H20,  AND  He 


by 


F.  Stuhl  and  H.  Nilci,  Scientific  Research  Staff 
Ford  Motor  Company,  Dearborn,  Michigan  48121 


Long  Abstract 

During  the  course  of  investigation  of  elementary  reactions 

. . 1 ° 

involving  OH  radicals,  we  have  flash  photolyzcd  (i  ^ 2 visec,  t > 1050  A) 

various  mixtures  of  CO  (0  to  0.2  torr) , H20  (0.05  to  0.25),  He  (20  to  170) 

in  the  presence  and  absence  of  02  (0  to  0.2).  A schematic  diagram  of  the 

pulsed  photolysis  apparatus  is  shown  in  Fig.  3.  The  kinetic  behavior  of 

OH  radicals  was  monitored  by  resonance  fluorescence  of  OH  (A2E+; 

2 1 

v*  = 0 -►  X it;  v"  = 0).  In  these  experiments,  the  initial  concentrations 
of  OH,  H,  and  0 were  estimated  to  range  from  2 to  30  x 1011  cm-3,  depending 
on  the  pressures  of  H20  and  02 . 
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Typical  results  are  shown  in  Fig.  4.  The  estimated  initial 
concentrations  of  OH  in  this  Figure  were  5 x 10**  for  the  runs  without  O2 
(Curves  A and  D)  and  1.5  x 10 12  cm-3  for  the  runs  with  O2  (Curves  B and  C) . 
The  diluent  He  was  kept  at  20  torr  in  these  runs.  Curve  A displays  the 
decay  of  OH  in  H2O  at  0.1  torr.  When  0.1  torr  O2  was  added  to  this  system 
(Curve  B) , the  lifetime  of  OH  was  shorter.  The  OH  decay  shown  in  Curve  D 
was  obtained  by  photolyzing  0.1  torr  H2O  in  the  presence  of  0.1  torr  CO 
but  in  the  absence  of  O2 . This  decay  can  be  explained  by  the  first-order 
removal  of  OH  in  the  reaction  OH  + CO  ->  CO2  + H.1  On  the  other  hand,  a 
significant  regeneration  of  OH  was  observed  when  0.1  torr  02  was  added  to 
the  above  mixture  of  0.1  torr  H2O  and  0.1  torr  CO  as  shown  by  Curve  C. 

I ? 

Ihis  behavior  could  be  similar  to  that  briefly  described  by  Greiner  in  the 
photolysis  of  a mixture  of  H202,  CO,  and  Ar. 

A detailed  analysis  of  these  data  and  additional  experimental 
work  are  under  way.  In  particular,  the  regeneration  of  OH  by  reactions 
involving  H02  such  as  H02  + O -*•  OH  + 02 ; H02  + H -*■  20H;  and  H02  + CO  -► 

CO2  + OH  are  being  examined. 

(1)  F.  Stuhl  and  H.  Niki,  "Pulsed  Vacuum-UV  Photochemical  Studies  of 
Reactions  of  OH  Using  Resonance  Fluorescent  Detection  Method," 

10th  Informal  Conf.  Photochemistry,  Stillwater,  Okla.,  May  1972, 
Abstract  E-4 . 

F.  Stuhl  and  H.  Niki,  "Pulsed  Vacuum-UV  Photochemical  Study  of 
Reactions  of  OH  with  H2,  D2 , and  CO  Using  a Resonance-Fluorescent 
Detection  Method,"  J.  Chem.  Phys.,  in  press. 

F.  Stuhl  and  H.  Niki,  "A  Flash  Photochemical  Study  of  the  Reaction 
OH  + NO  + M Using  Resonance  Fluorescent  Detection  of  OH,"  J.  Chem. 
Phys.,  in  press. 

(2)  N.  R.  Greiner,  J.  Chem.  Phys.  46,  2795  (1967). 
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OH  Concentration  (Arb.  Units) 


Fig.  3.  Schematic  diagram 
of  the  apparatus. 


Fig.  4.  OH  decays  in  the  H20-C0-He 
system  in  the  presence  and 
absence  of  63.  See  text. 
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CRITIQUE  OF  MEASUREMENT  TECHNIQUES 


During  the  first  afternoon,  there  was  considerable  discussion  of  the 
details  of  laboratory  techniques  . There  emerged  the  following  points  of  inter- 
est to  users  of  kinetics  data: 

1 . The  uncertainties  quoted  by  an  experimenter  usually 
are  based  on  the  scatter  in  his  own  data,  rather  than 
comparison  with  the  results  of  other  experimenters  . 

2.  In  some  cases,  the  evaluator  of  kinetic  data  is  also 
one  of  the  experimenters  whose  work  is  being  evaluated. 

3.  The  reactions  of  O atoms  are  fairly  well-measured 
compared  to  other  families  of  reactions  , and  the  results 
by  different  techniques  agree  reasonably  well  by  older 
standards . Yet  this  agreement  is  no  longer  good  enough 
for  some  of  the  central  reactions  in  ozone  aeronomy. 

4 . Any  reported  measurement  that  is  not  accompanied  by 

experimental  details  cannot  be  evaluated  and  should  not 
be  used.  The  most-quoted  example  is  the  Langley- 
McGrath  measurement  of  HO  +03"*  HO2  + O2,  which 
has  been  used  to  dismiss  odd-hydrogen  chemistry  from 
the  ozone  picture.  The  NBS  evaluators  feel  that  no  re- 
liable estimate  of  the  rate  of  this  reaction  can  be  made 
from  the  existing  data.  Consequently,  the  possibility 
remains  that  water  vapor  destroys  (or  creates)  ozone. 

E dels  on  pointed  out  that  the  modeler  can  help  the  kineticist  in  the 
design  and  interpretation  of  laboratory  experiments.  Some  kineticists  were 
receptive  to  this  suggestion,  feeling  that  free -radical  systems  are  getting  too 
complex  for  intuitive  reasoning.  But  other  kineticists  were  skeptical,  be- 
cause too  many  of  the  needed  input  rate  constants  are  uncertain;  at  this  stage, 
one  should  design  an  experiment  so  it  asks  only  one  question.  Gelinas 
assured  us  that  the  help  of  modelers  is  available  to  those  kineticists  who 
seek  it. 

REACTION  RATE  DATA  NEEDS 

H.  Johnston  stated  that  the  following  chemical  rates: 

HO  + NO2  + M "♦  HNO3  + M 
HNO3  + H + NO3 
HO  + NO2 

HNO3  + HO-*  H2O  + NO3 
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and  also  stratospheric  diffusion  coefficients,  are  needed  to  understand  the 
fate  of  odd  nitrogen  in  the  stratosphere.  Also,  the  temperature  dependence 
of  NO2  + O3  ■+  NO3  + O2  is  needed  to  establish  the  importance  of  the  NO3 
catalytic  cycle  for  destruction  of  ozone,  relative  to  the  NO2  cycle. 

H.  Schiff  emphasized  the  importance  of  measuring  the  wavelength 
dependence  of  the  quantum  yield  of  O('D)  in  ozone  photolysis  near  310  nm. 

The  atmospheric  yield  is  the  product  of  the  steeply-falling  quantum  yield  with 
the  steeply-rising  solar  intensity  in  this  region,  and  a lot  of  stratospheric 
chemistry  begins  where  O('D)  reacts  with  H2O  or  N2O.  The  relative  contribu- 
tion of  the  210  nm  region  to  O('D)  production  should  be  considered  by  modelers 
as  a function  of  altitude. 

There  is  a need  for  better  data  on  the  quantum  yield  for  photolysis  of 
N2O.  The  disagreement  by  a factor  of  2 among  various  experimenters  for  the 
rate  of  NO2  + O “*  NO  + O2  needs  to  be  resolved. 

Less  uncertain,  but  of  central  importance  to  the  ozone  problem,  is 
the  rate  of  O + O2  + M O3  + M . At  stratospheric  temperatures  , the  rate  ex- 
pression of  Huie  et  al  is  the  best  one  we  have. 

Further  refinement  of  the  rates  of  reaction  of  O('D)  with  H2O  and  NO2 
was  recommended. 

Cvetanovic,  representing  several  kineticists  , proposed  a systematic 
approach  to  the  collection  and  presentation  of  kinetic  data  . This  is  a matrix 
whose  row  and  column  headings  are  the  relevant  stratospheric  chemical  species. 
Each  square  contains  the  relevant  2-body  and  3-body  rate  expressions  and  in- 
formation on  their  reliability  (evaluated  results  of  several  experimental  studies, 
or  undocumented  guesses  by  experts)  , or  a statement  that  the  reaction  is  endo- 
thermic. One  kineticist  collects  and  assesses  all  reactions  of  a given  species. 
Because  of  the  huge  size  of  such  a matrix,  it  might  just  contain  references  to 
rate  expressions  listed  elsewhere. 

Johnston  accepted  this  format  for  inclusion  in  the  chapter  on  chemical 
reactions  of  the  CIAP  monograph  on  the  natural  stratosphere,  and  Garvin 
speaking  for  NBS , agreed  to  serve  as  the  assigner  and  collector. 

A similar  effort  for  UV  absorption  cross-sections  was  proposed,  per- 
haps updating  the  review  by  R.  D . Hudson  (Revs  . of  Geophys  . and  Space  Phys  . 
9^  305-406  (1971))  . With  regard  to  quantum  yields  , some  evaluations  are  in 
progress  at  NBS  , but  there  has  been  no  systematic  approach  to  all  molecules 
of  stratospheric  concern. 

Kummler  asked  how  CIAP  should  decide  which  minor  species  to  in- 
clude in  the  matrix.  The  most  important  species  are  those  that  influence 
either  ozone  or  the  major  infrared  radiators  . However,  the  web  of  uncertain 
chemical  interactions  implicates  many  others  , including  many  of  the  organic 
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species  containing  one  carbon  atom.  Also,  Ackerman  reminded  us  that 
species  whose  chemistry  is  understood  can  serve  as  tracers  to  estimate  eddy 
diffusion  coefficients;  those  so  used  already  include  CH4  and  N2O. 

UV  ABSORPTION  AND  SOLAR  FLUX 

Ackerman  discussed  the  Schumann -Runge  bands  of  O2,  and  gave  a 
cross-section  of  10“23  cm2  for  c>2  photodissociation  in  the  Herzberg  continuum. 
He  presented  ozone  spectra,  and  summarized  the  present  state  of  extraterres- 
trial flux  measurements.  The  6-bands  of  NO  are  shielded  by  the  Schumann- 
Runge  bands  of  O2  , whose  shape  is  temperature  dependent.  Thus,  a change  of 
30-40  K in  O2  temperature  can  change  the  NO  photodissociation  rate  by  one  or 
two  orders  of  magnitude  . 

Brasseur,  in  a theoretical  study  based  upon  the  above  data,  reported 
that  there  seems  to  be  no  net  production  of  NO  in  the  mesosphere  , and  there- 
fore no  downward  flux  of  NO  into  the  stratosphere. 

R.  D.  Hudson  found  predissociation  in  the  Schumann-Runge  bands,  in 
agreement  with  Ackerman . He  reported  that  the  unpublished  measurements  of 
Heath  show  no  variation  in  flux  at  250  nm  with  the  solar  cycle,  but  show  a 
variation  by  a factor  of  2.3  at  180  nm.  Hudson's  own  analysis  indicated  that 
a variation  by  a factor  of  3 was  possible  in  the  Schumann-Runge  bands  with 
solar  rotation,  and  that  this  would  affect  the  ozone  column  by  15%.  However, 
Ackerman  quoted  London's  recent  demonstration  that  there  is  no  correlation 
between  sunspots  and  atmospheric  ozone. 

J.  Park  described  his  calculation  of  the  mean  absorption  cross-section 
in  the  Schumann-Runge  bands  as  a function  of  altitude.  He  finds  little  varia- 
tion with  temperature;  the  analysis  is  appropriate  to  the  photodissociation  of  O2 
but  not  NO.  Hudson  and  Park  both  have  available,  for  use  by  others,  computer 
programs  for  the  computation  of  irradiance  and  O2  photodissociation  in  this 
region . 

Johnston  reported  the  absorption  spectrum  of  HNO3,  shown  in  the 
accompanying  figure.  The  spectrum  may  consist  of  two  or  more  electronic 
states,  with  dissociation  giving  H + NO3  in  the  shortwave  region  and  HO  + 

NO2  in  the  longwave  region.  Quantum  yield  studies  of  this  question  are  in 
progress  at  Berkeley.  Johnston  also  drew  the  spectra  of  NO2  (which  has  strong 
absorption  below  250  nm  besides  its  longwave  spectrum),  N2O5,  and  HN02* 

He  emphasized  the  importance  of  the  210  nm  window  to  stratospheric  chemistry. 
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Light  Absorption  Cross  Section  (cm  ) 


FIGURE  5.  ABSORPTION  SPECTRUM  OF  HNO3 
(Johnston  and  Graham,  to  be  published) 
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SOLAR  FLUX  MEASUREMENT  NEEDS 


Atmospheric  ultraviolet  flux  measurements  are  critical  to  the  success 
of  CIAP;  and  they  should  be  performed  simultaneously  with  the  measurement  of 
species  of  photochemical  interest.  However,  there  was  disagreement  on  the 
best  approach,  given  CIAP's  limited  time  and  resources  . It  was  asserted  that 
the  flux  in  each  narrow  wavelength  interval  can  be  computed  from  total  ozone 
column  measurements,  the  extraterrestrial  flux,  and  the  known  spectra  of  O2 
and  O3  . That  is  true  if  those  are  the  only  absorbing  species  and  their  tempera- 
ture and  pressure  dependence  is  insignificant.  This  should  be  verified  by  high- 
resolution  atmospheric  spectroscopy  near  210  and  310  nm.  Also,  good  measure- 
ments of  absolute  solar  flux  as  a function  of  altitude  are  needed. 

Many  of  the  above  desires  should  be  fulfilled  by  Heath's  planned 
aircraft -borne  measurements  of  solar  flux  at  1 nm  resolution  from  230  to  400  nm. 
But  in  addition,  the  Brewer  measurement  of  210  nm  flux  versus  altitude  should 
be  repeated  for  a variety  of  stratospheric  conditions  . 

Inexpensive  total  ozonesondes  should  be  flown  on  all  platforms  that 
are  measuring  other  species  . It  was  agreed  that  filters  with  10  nm  bandpass 
are  almost  useless  near  the  310  nm  ozone  absorption  edge,  because  of  its 
steepness . 

For  the  NO  predissociation  between  the  Schumann-Runge  bands  of  O2, 
spectral  measurements  of  better  than  0.01  nm  resolution  in  the  180  nm  region 
are  needed,  but  may  be  difficult  to  achieve  in  the  near  future. 


3-30 


IV.  SUMMARY  PROCEEDINGS  OF  THE 
ATMOSPHERIC  MODELING  PANEL 

(Thomas  D.  Taylor,  Chairman) 


INTRODUCTION 

In  April  1972  , the  first  CLAP  workshop  on  computational  modeling  of 
the  atmosphere  was  held  at  Asilomar,  Pacific  Grove,  California.  The  intent 
of  this  workshop  was  to  establish  the  state  of  the  art  in  computer  modeling  of 
atmospheric  flow  problems  which  relate  directly  to  the  CIAP  program.  The  pro- 
ceedings of  this  workshop  were  published  in  July  19  72  as  a Department  of 
Transportation  report  (TST-90-1) . 

As  a result  of  the  first  workshop,  the  CIAP  modeling  efforts  were 
divided  into  three  principal  parts: 

1 . Models  for  the  near  wake  flow  of  a supersonic  aircraft 
flying  in  the  stratosphere.  These  types  of  models  are 
limited  to  the  region  where  the  fluid  mechanics  is  domi- 
nated by  the  aircraft  itself  and  none  of  the  natural 
stratospheric  flow  phenomena  enter  the  model. 

2.  Transport  and  dispersion  models  which  predict  the 
spreading  of  the  aircraft  wake  in  the  stratosphere  for 
wake  dimensions  ranging  between  100  meters  to  1000  km. 

These  models  require  the  introduction  of  the  stratospheric 
flow  as  the  dominant  flow  mechanism. 

3.  Global  circulation  models  which  predict  the  changes  in 
climate  produced  by  the  dispersion  of  aircraft  emissions 
in  the  stratosphere. 
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The  first  workshop  arrived  at  the  following  conclusions: 

• Simplified,  approximate  methods  can  be  used  to  estimate 
high -altitude -aircraft  wake  dimensions  and  species 
concentrations  . 

• Two-  or  three-dimensional,  time -dependent  turbulent- 
flow  models  are  needed  to  treat  the  dispersion  and 
transport  region . The  models  should  include  provision 
for  nonequilibrium  chemistry.  It  should  be  determined  if 
radiation  is  important  on  this  scale  . 

9 Numerical  methods  appear  adequate  for  modeling  of  the 
wake  and  dispersion  and  transport  regions  . 

• Turbulent  flow  models  should  consider  the  effects  of 
turbulence  upon  reaction  rates  . Second-order  closure  and 
data  can  be  utilized  to  bound  the  magnitudes  of  eddy 
viscosity.  Application  of  second -order  closure  except 

to  estimate  effects  of  inhomogeneous  mixing  upon  chemical 
reaction  rates  presently  appears  impractical. 

m Presently  available  global  circulation  models  require 
major  refinements  and  extensions  to  be  applicable  to 
the  CIAP  problem. 

• It  is  unlikely  that  a three-dimensional,  global-circulation 
model  with  fully  interactive  chemistry  and  fluid  dynamics 
can  be  developed  by  the  19  74  CIAP  deadline. 

• The  suitability  of  various  simplified  models  for  treating 
global-scale  atmospheric  motions  and  physics  requires 
further  investigation . 

• Attention  must  be  devoted  to  development  of  a methodology 
for  coupling  models  of  the  wake  region,  the  dispersion 
and  transport  region,  and  the  global  circulation  in  order 

to  examine  effects  of  high-altitude  aircraft  operations  . 

9 Studies  are  needed  to  specify  accuracy  requirements  for 
various  atmospheric  models  applicable  to  CIAP  and 
to  compare  model  capabilities  with  required  accuracies. 

• Continual  and  expanded  interaction  between  fluid - 
dynamic  and  chemical  modelers  is  essential.  Furthermore, 
close  coordination  is  needed  between  CIAP  atmospheric 
modeling  investigators  and  other  CIAP  participants  (e.g.  , 
those  making  engine  emissions  measurements,  those  pro- 
jecting 1985-1990  high-altitude  aircraft  routes , those 
making  atmospheric  measurements  and  those  performing 
laboratory  chemistry  measurements) . 
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With  the  above  conclusions  as  guidelines,  the  atmospheric  modeling 
panel  was  charged  with  examining  the  various  CLAP  modeling  studies  to  deter- 
mine: (1)  if  the  current  efforts  will  provide  all  of  the  necessary  results  for 
CIAP;  (2)  how  the  modeling  efforts  will  interact  with  each  other;  and  (3)  what 
experimental  and  theoretical  information  is  required  for  input  conditions  to 
each  model. 

ATMOSPHERIC  MODELING  PANEL  PROCEEDINGS 
Near-Wake  Modeling 

The  panel  invited  the  various  investigators  to  discuss  their  efforts  in 
the  CIAP  program.  The  first  subject  discussed  was  near-wake  modeling. 

R.  Conti  from  Lockheed  Research  Laboratories  was  the  principal  speaker  in 
this  area.  Conti  discussed  the  Lockheed  model  for  the  flow  in  the  immediate 
vicinity  of  the  engine,  the  flow  in  the  vortex  interaction  region  and  the  effects 
of  buoyancy  and  shear  of  the  stratosphere.  The  details  of  the  model. are  avail- 
able in  a Lockheed  report  (LMSC-D309045)  and  therefore  will  not  be  presented 
here.  In  summary,  the  panel  members  concluded  that  the  Lockheed  fluid  flow 
model  for  the  aircraft -dominated  wake  region  was  adequate  for  CIAP  studies  . 

The  chemistry  model  in  this  region  was  only  briefly  discussed.  (It  was 
discussed  by  the  chemical  dynamics  panel.)  The  principal  limitation  in  the 
chemistry  results  seemed  to  be  the  uncertainty  in  the  input  data  for  the 
model.  The  results  of  CIAP's  engine  emissions  studies  should  reduce  these 
uncertainties  . 

Vortex -Interaction  Region 

In  the  vortex -interact ion  region  of  the  flow,  it  was  generally  agreed 
that  a detailed  flow  model  was  not  warranted  for  CIAP  since  in  practice  the 
chemistry  of  the  exhaust  products  is  essentially  frozen.  In  this  region, 
therefore,  only  a simple  model  which  accepts  the  species  from  the  very  near 
wake  and  convects  them  in  a frozen  state  to  the  dispersion  region  is  needed. 
Also,  due  to  the  short  time  of  exposure,  the  atmospheric  influences  of  radiation 
and  shear  can  be  neglected  in  the  near  wake . 

Dispersion  and  Transport  Region 

When  the  dispersion  and  transport  region  is  encountered  (at  times  of 
minutes  and  wake  dimensions  of  approximately  100  meters)  , flow  calculations 
become  dominated  by  atmospheric  shear  and  diffusion.  In  CIAP,  the  Lawrence 
Livermore  Laboratory  and  The  Aerospace  Corporation  are  developing  models  for 
this  region  of  the  flow.  The  work,  however,  had  not  progressed  to  a point 
where  detailed  presentations  were  possible.  J.  Walton  from  Livermore  did, 
however,  discuss  an  initial  study  on  "Scale  Dependent  Diffusion."  (A  con- 
densed version  of  Walton's  presentation  appears  later  in  this  chapter.) 

Walton  pointed  out  that  there  is  a lack  of  information  regarding  species 
diffusion  coefficients  in  the  stratosphere.  In  addition  to  Walton,  Y.  Pao  from 
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Flow  Research  and  G . Hilst  from  Aeronautical  Research  Associates  of 
Princeton  commented  on  the  dispersion  problem.  Hilst  pointed  out  that  the 
coupling  of  turbulence  with  the  chemical  reactions  can  significantly  affect 
the  prediction  of  reaction  rates  in  the  dispersion  and  transport  region.  The 
net  effect  is  a change  in  magnitude  of  the  laboratory -produced  reaction  rate 
constant.  Pao  pointed  out  that  the  magnitude  of  this  effect  can  be  estimated 
exactly  by  performing  detailed  calculations  of  a chemical  reacting  flow  using 
a method  developed  by  Dr.  S.  Orsag  of  MIT. 

Chemical  Modeling 

Following  the  discussion  of  the  intermediate  dispersion  and  transport 
region,  a presentation  was  made  by  J.  Chang  of  Livermore  in  which  the  subject 
of  computation  of  stratospheric  chemistry  was  considered.  Dr.  Chang  discussed 
the  numerical  integration  of  the  one-dimensional  species  transport  equations 
for  a restricted  chemistry  model.  From  his  study,  Dr,  Chang  concluded  that 
numerical  computation  of  stratospheric  species  chemistry  is  sensitive  to  the 
accuracy  of  the  numerical  fit  of  vertical  diffusion  coefficients  , In  addition, 
Chang  pointed  out  that  in  the  solution  of  the  stratospheric  chemistry  problems 
one  encounters  the  classical  "stiff"  equations  due  to  the  reactions’ approaching 
equilibrium.  As  a result,  a special  technique  must  be  employed  for  numerical 
integration  of  the  rate  equations  . 

During  Dr.  Chang's  presentation,  various  members  of  the  panel  had 
questions  regarding  the  proper  boundary  conditions  for  a chemistry  model  for 
the  stratosphere.  As  a result  of  an  open  discussion,  it  became  clear  that 
additional  experimental  information  is  required  in  order  to  reduce  the  uncer- 
tainty in  the  boundary  conditions.  J.  London  from  the  University  of  Colorado 
made  additional  comments  to  substantiate  this  conclusion. 

Global  Models 

The  discussions  of  global  models  by  the  panel  were  not  as  detailed 
as  the  other  models.  Descriptions  of  two  approaches  to  developing  global 
models  were  presented  by  M . MacCracken  from  Livermore  and  D . Cunnold 
from  MIT.  MacCracken  described  a simplified  "orange  slice"  global  flow  model 
for  the  stratosphere.  The  finite -difference  model  considers  variations  to  be 
minimal  along  latitude  lines  with  the  major  changes  occurring  in  the  vertical 
and  along  longitude  lines  . Cunnold  described  progress  in  developing  a spec- 
tral approach  to  modeling  the  global  circulation.  The  approach  is  to  approxi- 
mate variation  along  latitudes  and  longitudes  by  expansion  of  the  flow  variables 
in  spherical  harmonics.  The  coefficients  of  these  expansions  are  then  con- 
sidered to  be  functions  of  time  and  vertical  distance.  Equations  for  the 
coefficients  are  obtained  by  substituting  the  expansions  in  the  flow  equations 
and  performing  the  appropriate  integrations.  By  pursuing  this  approach,  the 
global  modeling  is  reduced  to  solving  a set  of  coupled  one -dimensional 
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unsteady  problems  . This  set  will  be  solved  by  a finite -difference  method. 
Cunnold's  presentation  will  be  published  as  a progress  report  to  the  CLAP 
Program  Office.  Due  to  the  limited  results  available  from  the  two  models  , 
the  panel  did  not  draw  any  firm  conclusions  regarding  the  usefulness  of  the 
approaches  . 

The  concluding  presentation  on  global  models  was  by  R.  Rapp  of  the 
Rand  Corporation.  He  discussed  Rand's  plans  to  perform  sensitivity  studies 
on  the  UCLA  global  circulation  model. 

CONCLUSIONS  AND  RECOMMENDATIONS  OF 
THE  ATMOSPHERIC  MODELING  PANEL 

After  the  global  modeling  presentations,  a summary  session  was 
held  in  which  the  panel  members  were  asked  by  the  chairman  to  review  the 
results  of  the  sessions  and  arrive  at  a set  of  recommendations  and  conclusions  . 

The  major  conclusions  of  the  modeling  panel  are: 

1.  The  near-wake  fluid  mechanics  calculation  of  Lockheed, 
is  adequate  for  the  CIAP  program.  The  chemistry  model 
for  this  region  of  the  flow  should  be  evaluated  carefully, 
however,  since  the  remaining  CLAP  modeling  studies 
utilize  the  output  of  these  predictions  . 

2.  The  details  of  the  fluid  mechanics  in  the  vortex-wake 
interaction  region  are  not  required  for  the  CIAP  study  since 
the  wake  chemistry  is  frozen  in  this  region. 

3.  Interaction  of  the  wake  with  the  stratospheric  radiation, 
shear  and  chemistry  can  be  neglected  for  the  near-wake. 

4 . The  effect  of  aerosols  on  the  near  wake  including  the 
vortex-wake  interaction  is  unknown. 

5.  Modeling  of  the  intermediate  dispersion  and  transport 
region  should  include  consideration  of  the  effects  of 
stratospheric  radiation,  chemistry,  and  shear. 

6.  The  effect  of  turbulence  on  the  chemical  reaction  rates 
should  be  evaluated  for  the  dispersion  and  transport 
regime  as  well  as  for  the  global-scale  regime. 

7.  Diffusion  coefficients  of  species  in  the  undisturbed  stra- 
tosphere appear  to  require  further  investigation  to 
support  transport  model  development. 

8 . The  effect  of  aerosols  in  the  dispersion  and  global 
models  requires  further  investigation. 
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9.  The  appropriate  boundary  conditions  for  species  transport 
models  in  the  stratosphere  are  uncertain. 

10.  The  radiation  and  chemical  kinetic  models  for  use  in  the  CIAP 
dispersion  and  transport'  regions  need  better  definition. 

11  . The  simplified  numerical  procedures  for  both  the  global  and 
wake  dispersion  calculations  are  still  not  fully  established. 

12.  Global  models  may  not  be  able  to  answer  CIAP  questions  by 

1974. 

13.  Chemistry  dispersion  models  should  be  uncoupled  from  the 
fluid  models  and  experimental  data  should  be  used  when  possible  to  provide 
shear,  diffusion  coefficients,  reaction  rates,  and  boundary  conditions  for 
CIAP  modeling. 

14.  Error  bounds  should  be  provided  by  each  modeler  regarding  the 
errors  in  the  input  data  as  well  as  the  model  itself. 

SCALE  DEPENDENT  DIFFUSION  1/ 

1.  Introduction 

G.K.  Batchelor  (1950),  in  studying  the  diffusion  of  a cloud 
about  its  center  of  mass,  locked  at  those  scales  of  motion  for 
which  Kolmogoroff  's  second  similarity  hypothesis  would  be  appli- 
cable. He  found  that  the  dispersion  of  a pair  of  particles, 
averaged  ewer  many  trials,  could  be  written  as  a simple  function 
of  time  for  limiting  cases.  At  an  early  stage  the  average  dis- 
persion was  quadratic  in  time,  while  later  it  became  cubic,  the 
transition  time  being  a function  of  initial  separation  oq  and  c, 
the  viscous  dissipation  rate.  In  this  note  it  will  be  shewn 
that  these  results  can  be  derived  from  a simple  diffusion  equation 
with  coefficients  proportional  to  the  4/3  power  of  the  cloud  size. 


— Material  presented  here  is  taken  from  a preprint  distributed  by  J.  Walton 
at  the  Atmospheric  Modeling  Panel.  The  full  paper  is  to  be  submitted  to 
the  Journal  of  Applied  Meteorology. 
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2.  The  Model 


Consider  a cloud  of  material  characterized  by  its  particle 
concentration  q(x,y,z,t)  with  q(x,y,z,o)  = a^(x,y,z).  Under  the 
action  of  turbulent  eddies  the  cloud  particles  will  be  assumed 
to  spread  subject  to  the  diffusion  equation 

H . Vt)  ^ . Ky(t,  0 . Kz(t,  0 


x,  y,  z being  coordinates  of  a frame  in  which  the  center  of  mass 
of  the  cloud  is  at  rest.  The  diffusion  coefficients  are  time  * 
dependent  through  their  relation  to  cloud  size  ox. 


Kx(t)  = cc 


1/3 


>4/3 

°x(t) 


(2) 


with  taken  to  be  the  standard  deviation.  The  constant  c is  of 
order  unity  arnd  may  contain  a factor  reflecting  the  geometry  of 
the  coordinate  system.  The  factor  e,  here  assumed  constant*  is 
usually  taken  to  be  the  atmospheric  dissipation  rate,  though  other 
parameters  have  been  proposed,  C.C.  Lin  (1960).  The  significance 
of  this  form  for  K lies  in  the  H/3  power  law  which  was  first  postu- 
lated by  Richardson  (1926)  on  the  basis  of  observations  available 
at  that  time  and  which  has  since  been  supported  by  dimensional 
arguments,  G.I.  Taylor  (1959).  It  is  tdnis  power  law  which  will  give 
rise  to  mean  square  cloud  dimensions  which  become  cubic  in  time  at 
late  time.  Since  the  relations  for  the  y and  z coordinates  are 
similar,  they  will  be  emitted  and  coordinate  subscripts  dropped 
where  possible. 

A general  solution  to  (1)  is 


q(x,y,z,t) 


(3) 


with  the  t's  given  by 


K(t')dt' 


(4 


For  a process  governed  by  (1)  the  variance  a obeys  the  equation 


gj-  = 2K(t) 


(s: 


Substitution  of  (2)  into  (5)  then  gives  the  equation 


do2/3  2 1/3 

-3F~’  1“ 


(6) 


with  solution 


a2  = (o2/3 
o 


2cc1/3t)3 


(7) 


= o2  ♦ 2a4/3ce1/3t  ♦ V/3c2c2/3t2  ♦ 
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where  a is  the  standard  deviation  at  t=o. 

© 

In  (7)  the  linear  term  dominates  on  the  interval 

1-1/32/3  t < 3 c-l/3q2/3  m<i  fran  its  coefficient  it  is  seen 
7c  o 7c  c 

that  observations  made  over*  short  time  periods  vri.ll  indicate  Fickian 
1/3  ii/3 

diffusion  with  K » ce  oq  . Cn  the  interred 

c“^3o2/3  < t < e^c2'3  o2  is  quadratic,  while  there- 

after growth  is  dominated  by  the  cubic  tern.  The  coefficients  of 

2 3 . . 

t and  t ®id  the  time  of  transition  are  in  agreement  with  those 

of  Batchelor  (1950). 


The  reader  will  be  aware  of  the  arguments  for  treating  large 
scale  atmospheric  motions  as  characteristic  of  turbulence  in  two 
rather  than  three  dimensions.  The  added  constraint  of  two  dimen- 
sional turbulence  gives  rise  to  a diffusion  coefficient  prcpor- 
2 4/3 

ticnal  to  o rather  than  o , which  upon  substitution  into  (5) 
gives  cloud  growth  exponential  in  time.  This  problem  has  been  dis- 
cussed in  a reoent  paper  by  Jung-Tai  Lin  (1972)  and  will  not  be 
explored  further  here. 
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3.  Dctensi.cn  to  very  late  time 

Vhile  no  prediction  car  be  made  cn  the  basis  of  the  above 
equations,  it  is  knam  (Ch.  3,  F.  Pasquill  (1962))  that  after  sore 
tine  tp  a function  of  the  properties  of  the  flew,  growth  character- 
istic of  a constant  diffusion  coefficient  will  again  be  observed. 
After  t^  Lq.  (3)  for  the  concentration  q is  still  correct  while 
(4)  simplifies  to 


o 


Cloud  size  is  obtained  by  substitution  of  constant  in  (5)  giving 


a2  = o*  ♦ 2K1(t-t1) 


(7*1 


where  is  the  cloud  size  at  t * t^. 

4.  Comparison  'with  Observation 

In  Fig.  1,  standard  deviations  obtained  from  (7)  and  (?') 
in  one  dimension  are  carpared  with  an  observed  growth  curve  (the 
small  circles  in  this  figure)  presented  by  Randerscn  (1972).  In 
his  work  Randerscn  analyzed  plain-view  diagrams  obtained  by 
Laurence  Livermore  Laboratory  from  observations  of  the  debris 
cloud  resulting  from  a nuclear  cratering  experiment.  For  the 

purposes  of  comparison,  the  dissipation  rate  e used  in  (7)  was 
2 3 

2.6  an  /sec  , from  values  at  700  irb  and  500  irb  calculated  by 
Ellsaesser  (1969).  The  constant  c was  1.0  and  the  initial  size 
oq  was  taken  to  be  0.65  km.  Curve  I exhibits  growth  for  which 
no  upper  limit  on  K is  imposed.  After  1-1/2  hrs  this  growth  is 

dominated  by  the  cubic  term  in  (7).  For  curve  II,  K had  an  upper 
. . 10  2 

limit  of  2x10  cm  /sec  corresponding  to  values  of  K found  through 

the  study  of  large  scale  motions  in  the  atmcsphere,  f*urgatrcyd 

(1969).  Curve  III  was  obtained  using  Randerscn' s asymptotic 
8 2 

value  of  K,  7x10  cn  /sec,  hence  "die  very  pood  agreement  in  this 
range.  Curves  II  and  III  illustrate  the  wide  range  of  values  of 
K found  from  different  types  of  observations  and  the  correspond- 
ing problem  of  describing  cloud  dispersion  in  this  regime.  The 
most  apparent  difference  between  the  theoretical  and  observational 
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curves  is  in  the  irore  shallow  slope  of  the  latter.  It  is  felt 
that  this  is  due  to  removal  mechanisms  which  are  not  present  in 
the  simple  model  described  here. 
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Fig.  6.  Comparison  of  computed  standard  deviation,  o(km)  with 
observation  (small  circles),  Randerscn  (1972):  I.  Ho  upper  limit 

cn  diffusion  coefficient;  II.  Maximum  diffusion  coefficient 
Kx  = 2x10 10  cm2 /sec;  III.  = 7x10®  an2 /sec. 
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